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Abstract 

A direct measure of hydrogen bonding in water under conditions ranging from the normal state 
to the supercritical regime is derived from the Compton scattering of inelastically-scattered X- 
rays. First, we show that a measure of the number of electrons rig involved in hydrogen bonding 
at varying thermodynamic conditions can be directly obtained from Compton profile differences. 
Then, we use first-principles simulations to provide a connection between Ue and the number of 
hydrogen bonds nuB- Our study shows that over the broad range studied the relationship between 
ne and uhb is linear, allowing for a direct experimental measure of bonding and coordination in 
water. In particular, the transition to supercritical state is characterized by a sharp increase in the 
number of water monomers, but also displays a significant number of residual dimers and trimers. 
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The complexity of the phase diagram of water and its relevance to many physical, chemical 
and biological processes is to a large extent related to the presence of hydrogen bonds ^ . 
Over the years, different methods have been used to investigate the structure and dynamics 
of water [2]]; nevertheless, many fascinating questions still remain. Hydrogen bonding and 
coordination in disordered phases, most notably in the liquid, are the subject of intense study 
and the average number of hydrogen bonds {uhb) in water has recently become a subject of 
controversy |3|, X-ray absorption spectroscopy (XAS) [3] in water and ice indicated that 
the liquid contains significantly more broken hydrogen bonds than previously thought, but 
these conclusions have been questioned [4,5]. The degree of persistence of hydrogen bonding 
in the region beyond the critical point is also an interesting question Gj. Several works report 
on residual hydrogen bonding in supercritical water [2I, The number of hydrogen bonds 
per molecule nns is a concept of some ambiguity since it depends on geometric definitions. 
The connectivity of the hydrogen bonded network has also been recently studied in terms 
of electronic signatures {s] highlighting the presence of a fast electronic bond dynamics and 
reorganization. 

In this work, we highlight how Compton scattering - i.e. inelastic scattering of photons by 
electrons jsl - can actually measure the number of electrons involved in hydrogen bonding, 
and we use this result to characterize structure and bonding in water under a variety of 
thermodynamical conditions. The Compton profile (CP) is given by 

J(p,e) = j n(p')5(p'-e-p)dp', (1) 

where e is the unit vector along the scattering direction and n(p') is the electron momentum 
density; here we only consider the spherical average J{p) over e. The unique sensitivity of 
CPs to valence electrons and to chemical bonding (via ^(p')) is at variance with diffraction 
techniques, which are principally sensitive to ionic positions. Earlier work usin g C ompton 
scattering has in fact revealed the quantum nature of the hydrogen bond in ice 10|, and it 
las also been s ugg ested that CP can be a measure of the amount of H-bonding in water 
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17l |. The Compton profile reflects the distribution of electronic mo- 



mentum in the ground state of the measured system. Changes in bonding induce changes 
in this distribution and by tracking these changes both in experiment and calculations for 
water we establish two central points. First, we show that we can measure directly the 
number of electrons associated with change in bonding induced by varying thermodynamic 



conditions. Then, we demonstrate a very robust linear relationship between uhb and Ue, 
the number of electrons per molecule involved in hydrogen bonding, from ambient to super- 
critical conditions. Therefore, the measured rig becomes a direct and unambiguous measure 
of Uhb- 

As shown in FiglH our measurements cover the temperature range along the isobar P=300 
bar, and include the following state points: (a) T = 30 °C (density = 1.01 g/cm^), liquid 
water reference; (b) T = 200 °C (density = 0.89 g/cm^), T = 300 °C (density = 0.76 g/cm^) 
and T = 350 °C (density = 0.67 g/cm^), all showing precursor effects of hydrogen-bond 
breaking and change of local tetrahedral structure; (c) T = 397 °C (density = 0.40 g/cm^), 
critical point neighborhood; (d) T = 416 °C (density = 0.21 g/cm^), supercritical conditions. 

Our measurements were taken at the European Synchrotron Radiation Facility in Greno- 
ble at beamline ID15B using the 56 keV experimental set-up and a specially designed high- 
pressure cell [l8|. The scattering angle was set at 144.4° and the back-scattered photons 
for the CP were collected by a Ge multi-element detector. After subtracting the very weak 
background (5xl0~'^% at the Compton peak), the raw data were coverted from wavelength 
scale to momentum scale, and multiple scattering contributions were subtracted from the 
measured profile in order to obtain the total profile. The multiple scattering contributions 
were calculated using Monte-Carlo simulations with each dataset treated separately in or- 
der to account the change in density. The ratio of multiple scattering to the total ranges 
between 0.9% at 30 "C to 5% at 416 °C. Our statistics ranges between 2.1x10^ counts at 
30 °C to 1.3x10^ counts at 416 °C at the Compton peak. The energy-dependent resolution 
function was obtained from the full-width at half maximum of the elastic peak (0.51 atomic 
units of momentum). Now, since the CP is a projection of the electronic momentum density, 
its integral yields the total number of electrons and provides a convenient normalization. 
In this work, we focus our attention on CP differences, which are free from background 
contributions as long as the sample environment is unchanged. We limit ourselves to the 
stated thermodynamic range because measuring the CP of monomers (coordination 0) is 
not possible with the required high statistics and measuring ice Ih (coordination 4) would 
require a different sample environment. We thus present our data as the difference AJ(p) 
between a reference state CP (taken as liquid water at 30 °C and 300 bar) minus the CPs 
at higher temperatures. The acquisition times are of several hours. 

Our Compton scattering calculations rely on first-principles methods without any struc- 
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FIG. 1: Experimental measurements were taken along the isobaric curve for P=300 bar. 

tural information from experimental data (except for the density) [19]. To compare ex- 
)eriments with theory, we use extensive Car-Parrinello molecular-dynamics simulations 
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22], using 32 water molecules in periodic-boundary conditions. The subtle tech- 



nical issues arisin g in the simulations of liquid water have been discussed extensively in the 



literature 
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26| - in particular it has been suggested [27|] that a closer agreement 



with experiments is obtained by constraining the internal degrees of freedom (i.e. keep- 
ing each water molecule rigid), given that at the conditions considered those vibrational 
degrees of freedom are frozen out and consequently in their zero-point motion state. We 
thus performed first-principles molecular-dynamics simulations at ambient, subcritical and 
supercritical regimes at 77, 200, 300, 350, 397 and 430 °C [28| with rigid water molecules. 
To determine the average water molecular structure to be used in the simulations men- 
tioned above and to understand the effects of the change in intramolecular geometries on 
CP differences at various conditions, we also performed fully unconstrained Car-Parrinello 
molecular dynamics simulations [29\. These unconstrained simulations will also be used 
later to further confirm the linear relationship of rie with uhb- We further choose for the 
rigid-water molecular structure the average intramolecular 0-H distance {toh) and H-0- 
H angle (Zhoh) to be 0.995 A and 105.0° respectively, which are obtained from the 300 
°C simulation. We note as expected in our unconstrained simulations, a small change in 
average structure over the range of temperatures considered, with toh and ^hoh varying 
between 104.6 - 106.1° and 0.990 - 0.998 A respectively, when the temperature decreases 
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from 430 to 127 °C. These changes have, in principle, an effect on the CPs, but it will be 
shown that this effect is reduced to that of a scale factor: the linear relation alluded to 
above remains robustly verified. All our simulations were done in the NVT ensemble, with 
the densities fixed at the experimental values at 300 bar [30], using ultrasoft pseudopoten- 
tials with a plane-wave kinetic energy cutoff of 25 Ry for the wavefunctions and 200 Ry 
for the charge density. The wavefunction fictitious mass is chosen to be 1100 a.u. and a 
time step of 10 a.u. is used to integrate the electronic and ionic equations of motions. It 
is important to stress that these choices of fictitious mass and time step are the same as 
those used in Ref. and are accurate for this system and wavefunction cutoff - i.e. they 
provide a structure that is not affected by integration errors. We note that fully-flexible 
water molecules described with norm-conserving pseudopotentials would require smaller fic- 
titious masses (~340 a.u.) to provide accurate results - larger values will integrate less 
accurately the equations of motion, and lead to less structured radial distribution functions 
as compared to tightly-controlled simulations (such choices improve rather artificially the 
agreement with experiments, and lower the theoretical freezing point of water) 23 1. 

The theoretical reference state is chosen to be that at 77 °C: a lower temperature would 

n 

see the onset of a freezing transition [27|. The reason for the shift of the liquid-solid phase 



25l . l27l | ; increasing the 



boundary have been discussed extensively in the literature [23|, |24J . 
temperature for the only case (30 °C) close to the liquid-solid phase transition avoids the 
pitfall of remaining trapped in a glassy configuration. As a validation of the theoretical 
temperature scale, we found good agreement between the predicted and the measured self- 
diffusion coefficients. The predicted (and, in parenthesis, measured [3ll, l32|) self-diffusion 
coefficients obtained for water at 77 (30 in the experiment), 200, 300, 350, 397 and 430 °C 
are 1.5 (2.6), 9.5 (18.9), 27.0 (37.4), 35.3 (45.3), 75.9 (80.3) and 163.3 (174.4) xlO-^cmVs, 
respectively. We therefore adopt hquid (rigid) water at 77 °C, as a theoretical counterpart 
for the experimental temperature of 30 °C since the lowest temperature for which dynamics 

of water remains in the liquid phase is about 50 °C. ^ ^ 

To compute J{p,e), we used the approach by Romero et al. [33], based on maximally- 
localized Wannier functions j^]- Each water simulation lasted about 17 ps after thermaliza- 
tion, and at each temperature the CPs are calculated from the averages of 10 uncorrelated 
configurations, equally-spaced by 1.7 ps and then spherically averaged. The CPs are also 



convoluted with a Gaussian that matches the 0.51 a.u. experimental resolution. 
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FIG. 2: (color online), (a) Experimental CP difference AJ(p) (points, in electrons per atomic 
unit of momentum) with respect to the experimental reference state as a function of temperature, 
and theoretical CP difference A J(p) (solid lines, in electrons per atomic unit of momentum) with 
respect to the theoretical reference state as a function of temperature, and scaled by 0.73. The 
temperature of reference state is 30 °C for the experiments, 77 °C for the theory (at a density 0.97 
g/cm^). For the purple experimental points the temperature is 416 °C while for the purple theory 
curve T = 430 °C. (b) The number of electrons rie involved in hydrogen bonding (i.e. integrated 
absolute profile difference) as a function of temperature (simulations have been scaled by 0.73). 

We first show the experimental and theoretical A J(p) in Fig. [2]^a). The overall agreement 
between experiment and theory is excellent, provided the theoretical curves are rescaled to 
compensate for a systematic overestimation of the experimental values. The reasons for thi s 
overestimation can be manifold; they have been already reported in the literature [lO 
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33|, 



but will be discussed in detail later in this paper. What should be stressed however is that 
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a unique rescaling factor (here, 0.73) brings all theoretical curves in quantitative agreement 
with experiment. 

We now propose an electronic measure of hydrogen-bonding defined by the number of 
electrons 

n, = l r \AJ{p)\dp, (2) 



2„ 

where, as mentioned before, AJ(p) is the difference between the CPs at the conditions 
studied and the reference profile for normal liquid water. The area rig (in units of electrons) 
between a difference curve and the horizontal axis in Fig. [2] measures the amount of electrons 
whose wavefunctions change in going from the reference state to the one being measured. The 
quantity Ue tracks the progressive breaking of hydrogen bonds as the temperature is raised 
from ambient conditions. Both the experimental measurements for Ue and the theoretical 
predictions (after rescaling) are plotted in Fig. EJ^b) as a function of temperature. All the 
nf,{T) points line up, reflecting the agreement between the experimental and theoretical 
AJ{p) in Fig. [2](a). Interestingly, ne(T) in Fig. [2](b) does not saturate, indicating that even 
at the highest temperatures considered significant hydrogen bonding remains. Moreover, 
the almost linear behavior of np(T) in the range from 200 to 430 °C mirrors a similar trend 

n 

for uhb, as shown in Ref. [2] (Fig. 13). 

These considerations link rig to the variation in hydrogen bonding as the temperature 
of the system is changed. In the following, we make a direct connection between ng and 
^HB by studying a model cluster where these quantities are unambiguous: a five- molecule 
model with one molecule at the center and the other four surrounding it in a tetrahedral 
arrangement. The 0-0 distance between two water molecules in the cluster geometry is fixed 
at 2.77 A, which is the first peak height position of the theoretical 0-0 radial distribution 
function of liquid water at ambient conditions within our rigid water simulations. It has been 



shown 16j that the CP differences are strongly dependent on the 0-0 distance. An increase 
in the 0-0 distance will lead to a decrease in ng and this can affect the scaling factor, as 
we shall discuss later in the paper. In the present cluster model, a single hydrogen bond 
can be formed or broken simpl y by removing, one at a time, the four outer molecules. This 
model was introduced in Ref. [iJ], but here, using maximally-localized Wannier functions 
3J], we can resolve the CP unique to the central molecule, as the hydrogen bonds with its 
neighbors are broken, one at a time. The number uhb of hydrogen bonds formed by the 
central molecule is unambiguously defined, and the corresponding CP has been calculated 
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FIG. 3: (color online). The number of electrons rig involved in hydrogen bonding (scaled by 0.73) 
as a function of hydrogen bonds uhb- The zero on the vertical scale corresponds to the reference 
state at ambient conditions. Squares: water clusters (1 to 5 molecules); circles: rigid-water bulk 
simulations. 

by averaging over all the possible clusters with a given uhb- can thus be plotted as a 
function of uhb- The results are shown in Fig. [3] (squares), where rig is again rescaled by 
a factor of 0.73, to facilitate a later connection with experiments. A linear dependence for 
the cluster results is seen and the slope provides an indication of the number of electrons 



involved in one hydrogen bond 35 1 



Finally, our first-principles molecular dynamics results are also plotted on Fig. [3l using 



for Uhb the structural definition 



361] based on bond-length and bond-angle criteria. The 



{nejUHB) points collapse on the same straight line, confirming a universal relationship be- 



tween Ue and Uhb, and demonstrating that the criteria [36|] are meaningful and consistent 
with the natural integer uhb count of the cluster model. Remarkably, the linear relationship 
between Ue and uhb persists under a wide variety of conditions, ranging from the normal 
to the supercritical regime and covering all degrees of hydrogen bonding. Moreover, the 
slope in Fig. [3], as mentioned earlier, provides the number of electrons involved in a hy- 
drogen bond: (13 ± 0.6)xl0~^ electrons are displaced for every bond that is broken (or 
formed). The bonding and coordination picture that emerges from these results is also in 
close agreement with the established understanding of coordination in liquid water. Our 
results are summarized in Table [H where we can see that in going from the normal state 
into the supercritical regime, uhb decreases from 3.40 (compared to the value 3.58 from 
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FIG. 4: (color online). The (unsealed) number of electrons rig involved in hydrogen bonding as 
a function of hydrogen bonds uhb for flexible water simulations at 200, 300, 403 and 430 °C. 
Note that the reference state in these calculations is 127 °C due to freezing at lower temperatures 



24 



Ref. 6{]) to 0.34. Moreover, except at ambient conditions, the theoretical uhb^^ at different 
temperatures agree reasonably well with the NMR studies 2|], which also suggests that the 
simulations and experiments produce similar structures at the same temperatures above the 
theoretical melting point. 

We now tackle in detail the origin of the rescaling factor. In principle, it is possible to 
exactly reproduce the amplitudes of the experimental CP differences by adjusting the in- 
tramolecular 0-H distance and the H-bond geometry; this was detailed in Refs. 16| and 17|. 
Our approach uses instead geometries that are obtained directly from first-principles simula- 
tions and identifies a-posteriori the rescaling factor that aligns experiments and simulations. 
Note also that the proton zero-point motion is not accounted in any these simulations and 
may have a similar impact as discussed by Isaacs et al. flO] for the CP anisotropy of ice 
(the phenomenological rescaling used in those calculations was 0.6). For this same system. 



Romero et al. 



331 ] have shown that finite temperature effects can be ruled out as the reason 



for rescaling and Ragot et al. 121] demonstrated that correlations effects are also not crucial. 
It is beyond the scope of this work to determine whether the proton zero-point motion or the 
difference in geometries due to the approximate exchange-correlation functionals dominate 
the discrepancies between experiments and theory. Instead, we show that the linear relation 
between rte and uhb holds even when intramolecular distances are allowed to change with 
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FIG. 5: (color online), (a) Snapshots from the simulations at 200 °C, 300 °C and 430 °C. Each 
panel shows the 32 inequivalent molecules in the periodic cubic supercell (note that the volumes at 
these conditions are very different, and the figures presented here have been rescaled accordingly). 
The color codes the connectivity of each water molecule (see text). Blue: water molecules in 
clusters with 6 or more molecules; green: water molecules in clusters with 2 to 5 molecules; red: 
monomers, (b) Number of water molecules riw that belong to a cluster of size rid, at different 
temperatures. 

temperature in fully unconstrained simulations. Even though these simulations are closer 
to the experiment as far as scaling is concerned we prefer to use the rescaled rigid water 
calculations for the fundamental reasons given above. Fig. [3] and Fig. H] clearly show the 
remarkable robustness of our linear relation. Finally our fully unconstrained simulations 
clearly show that Compton profiles may be affected by changes in internal molecular geom- 
etry, but these changes are caused by hydrogen bonding and cannot be treated in a separate 
way. 
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We now examine the connectivity of the hydrogen-bond network as we move from the 
normal to the supercritical regime. For each molecule, we can determine the number of 
neighbors on a connected hydrogen-bond network. Fig. O^a) shows three snapshots taken 
from our rigid water simulations. We see that at 200 °C most molecules (in blue) are 
in clusters containing 6 or more molecules, and very few monomers (in red) exist. Most 
remarkably, no intermediate cluster size exist - i.e. a molecule is either fully connected to 
all other, or fully disconnected, probably reminiscent of the very high surface tension of 
water. At 300 °C, there is a very wide distribution of cluster sizes: monomers and small 
clusters containing 2 to 5 molecules (in green). Finally, at 430 °C monomers dominate, but 
a few clusters of less than 6 molecules persist. To quantify these observations, we show in 
Fig. MJo) the distribution of cluster sizes at a given simulation - molecules belong to a 
cluster (of size nd) if they are hydrogen bonded together. The uppermost panel of Fig. Mjo) 
shows that at normal conditions almost all water molecules are connected to each other. 
The percolating hydrogen-bond network breaks somewhere between 300 °C and 397 °C, as 
illustrated in Fig. [5](b) by the diffuse and bimodal distribution at 300 °C and the peaked, 
monotonically decreasing distribution at 397 °C. According to Table [H this is consistent 
with the well-known percolation threshold uhb = 1-53 — 1.55 |37^. Finally, the lowest panel 
in Fig [5](b) shows a distribution at 430 °C with a dominating amount of monomers but also 
a substantial amount of dimers and trimers, indicating that the amount of hydrogen bonds 
is still non-negligible in supercritical conditions. The experiment yields nHB=0.74±0.16 at 
416 °C, while the simulation predicts nHB=0.34: at 430 °C. These values are consistent with 
Fig. 13 in the review article by Kalinichev 2|. 

In conclusion, we have shown that Compton scattering provides an absolute measure 
of the number Ug of electrons per water molecule involved in hydrogen bonds. We have 
measured the dependence of function of temperature and compared with state-of- 

the-art first-principles molecular dynamics simulations. This comparison leads to a linear 
relation between Ug and the average number uhb of hydrogen bonds per water molecule. The 
linear relation is very robust and holds over all possible coordination numbers. This approach 
should prove valuable in the future to address issues that arise in the characterization of 
hydrogen-bonded networks. 

We acknowledge useful discussions with G. Loupias and V. Honkimaki, and thank B. 
Wood for his help and suggestions in calculating the connectivity of a bonding network. We 
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TABLE I: Experimental and theoretical (scaled by 0.73) rig, and theoretical uhb at different state 
points, uhb is calculated from the molecular dynamics trajectories using the structural criterion 
of Ref. Q]. The first column shows the experimental temperature and pressure at each state point. 



Temperature 


rip flO-2) 


rip flO-2) 


77. rr R 


and pressure 


(experiment) 


(theory) 


(theory) 


200 °C, 300 bar 


1.27 ± 0.19 


1.28 


2.46 


300 °C, 300 bar 


2.03 ± 0.19 


1.96 


1.71 


350 °C, 300 bar 


2.58 ± 0.20 


2.39 


1.38 


397 °C, 300 bar 


3.07 ± 0.21 


3.09 


0.73 


416 °C, 300 bar 


3.30 ± 0.21 






430 °C, 300 bar 




3.60 


0.34 
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